Tails of Localized Density of States of Two-dimensional Dirac Fermions 
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. The density of states of Dirac fermions with a random mass on a two-dimensional lattice is 

considered. We give the explicit asymptotic form of the single-electron density of states as a function 
of both energy and (average) Dirac mass, in the regime where all states are localized. We make 
' ■ use of a weak-disorder expansion in the parameter g/m^, where g is the strength of disorder and 

(~| ' m the average Dirac mass for the case in which the evaluation of the (supersymmetric) integrals 

, corresponds to non-uniform solutions of the saddle point equation. The resulting density of states 

^ ' has tails which deviate from the typical pure Gaussian form by an analytic prefactor. 
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I . I. INTRODUCTION 

a .... . 

Q ' Dirac fermions in two dimensions play a prominent role in the description of some notable two-dimensional (2D) 
^ O j condensed matter systems. Examples are the plateau transition in the Integer Quantum Hall Effect (IQHE) D- 
wave superconductivity in CuO-planes [|| and quasiparticles in the resonant valence bond state of the two-dimensional 
[ Heisenberg model Q]. Notwithstanding, of course, the 2D Ising model in the framework of a Grassmann field repre- 
I> . sentation (for a review, see |||] ) . The Dirac mass creates a gap between the hole and the particle band and, therefore, 
OO ' gives an opportunity for the description of an insulating fermionic system when the Fermi energy lies inside this gap. 
The tuning of the mass allows us to go through a transition to a metallic state when the mass (i.e. the gap) vanishes. 
This behavior is relevant for the description of, e.g., the plateau transition in the IQHE. However, at zero energy 
I (measured with respect to the Fermi level) the density of states (DOS) is always zero, as it turns out, for the pure 
• Dirac fermions, even for the case of a vanishing mass. 
^\ In practice, in a real system there are always relevant impurities such that the pure Dirac fermions model is 

insufficient as a description. In order to take these impurities into accout we can simply add randomness to the Dirac 
fermions' Hamiltonian. This has severe consequences on the properties of the Dirac fermions; in particular, as it turns 
^ ■ out, randomness is a non-perturbative effect Randomness can lead to new (localized) states even in the gap of 
I the pure system due to the formation of Lifshitz tails. This phenomenon will be studied in the present paper for the 
' ^ , limit case of a large mass. We stress that, although there are different types of randomness which can be in principle 
' included in the Dirac model, we will concentrate here on the physically relevant case of a random Dirac mass only. 
Q ] Moreover, we will only investigate here the form of the average DOS for the localised states. The average DOS has 
• • ■ been evaluated via various other approaches. Fradkin, for example, computed this quantity for Dirac fermions with a 
. 5^ ' random energy 0| by means of the coherent potential approximation (CPA) and found a filling of the pure system's 
I gap. In the case of a Lorentzian mass distribution, e.g. P{Mr) = (r^ -I- {Mr — m)^)~^r/7r, the average DOS can be 
' even calculated exactly At energy E = it has, in contrast to the vanishing DOS of the pure Dirac fermions's 
5^ , case, a non-zero value 

p{m,E = 0) = p{rn) = rlog[l + A/(r^ + rn'^)]. (1.1) 

The lattice constant 1/ A shows up as a cut-off for an infrared divergency. The power law of the tails at large mass 
m is an artefact of the Lorentzian distribution, which has strong tails itself. For more realistic distributions, like 
the Gaussian P{Mr) = exp(— (M^ — m)'^ /2g) / y/¥g, we expect much weaker tails; for this case the renormalization 
group (RG) calculation gives a vanishing DOS. This is a consequence of the fact that a random mass is an infrared 
marginally-irrelevant perturbation which drives the system always back to free Dirac fermions' theory as its fixed 
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point ||9|,|3|,|l0|,y_l| . This somewhat bizarre result must be traced back to the intrinsic perturbative nature of the RG 
calculation: the /3-function is indeed always evaluated by means of a perturbative expansion of the disorder strength 
g. Although the solution of the RG equations may pick up some non-perturbative feature for small g, since these 
equations make statements valid to all orders, there is nevertheless no guarantee that all non-perturbative effects are 
duly taken into account. There are in fact rigorous estimates of the DOS at m = which give a deviation from the 
vanishing DOS of the free Dirac fcrmions' fixed point 

p{0) > cie-^^/^ (1.2) 

where ci > and C2 > are constants independent of g. This result is similar to the one from the (uniform) 
saddle-point (SP) approximation (large N limit), where p(0) ~ e^^^^ /g was found 

For the Dirac fermions, a non-perturbative approach based on the SP approximation of a suitable functional- integral 
representation has given, again for the extended states, a finite bandwidth in the DOS having a characteristic semi- 
circular form I p^ . Perturbative calculations (leading, however, to different and sometime unphysical results) have 
also concentrated on the extended states. Therefore, hitherto no specific information on the features of the DOS 
for the localized states has been obtained from the Dirac fermion approach and it certainly of some great interest to 
characterize the localized states in view of a possible resolution of the still open problem of the localization length 
near the Integer Quantum Hall Transition (IQHT). 

In this paper, we produce such a calculation for the Dirac fermions model characterized by a Gaussian random 
mass distribution. In order to make the calculation feasible, we consider the behavior of the non-uniform SP solution 
in the small g/rri^ regime. The behavior for small g/m? can be in fact explained in terms of an SP approach to the 
supersymmetric functional integral for the averaged DOS [|3|jlJ]. Instead of doing the SP integration, we choose to 
expand the functional integral directly in powers of 1/m. 



II. THE MODEL 



In this Note we carry out an investigation on the form of the localized DOS at zero energy, which leads us to 
characterizing the tails of this quantity in the asymptotic limit of a large average mass (or, as it turns out, in the 
limit of weak disorder in the single-particle random potential). Our starting point is the Dirac fermions' Hamiltonian 
for the independent quasiparticles 

Hd = iVio-i -I- iV2CT2 -I- M(T3 = iV • (7 -I- Af(73, (2.1) 

where the energy is measured in units of the hopping parameter of some original lattice model, Vj is the lattice 
differential operator in the j-direction and {cj} are the Pauli matrices. This Hamiltonian, with a random mass term 
Mr and no random vector potential, is a reasonable starting point for the description of a number of interesting 
physical systems, as discussed in the Introduction. The average local DOS is obtained as usual from the averaged 
one-particle Green function G{E — ioj) = {Hd — E + iuj)~^ as 

pim,E) = - lim -Im(TT2GrriM,E -iuj)), (2.2) 

where Tr2 stands for the trace over the 2x2 matrix structure. We have characterized the random Dirac mass 
effectively by its average value (Mr) = m and correlation function {Mr Mr') — gSry- The representation for the 
Green's function in terms of both commuting and anticommuting functional integrals is well-known [ p^ and we arrive 
at a collective fields representation which for the Green's function with r' — r gives 

{Grr{M,E -iLj)) ^- I TQrTe-^VPVQVeVQ. (2.3) 
9 J 

T is here the diagonal matrix (1,«), i.e = a^. The mixed commuting-anticommuting (supersymmetric) effective 
action has the form (setting, for convenience, E — Q) 

S = i y (TraQ^ + Tr2P2 + 2Tr2e,e,) + 
9 

+ lndet[(i?o + it^o-o - 2tQt){Hq + iojcro + iirPry^] 

+ lndet[l - AtQt{H'^ + iojaa + 2iTPTy^TQT{Ho + iojcra - 2TQr)"i]. (2.4) 
Here, Hq is the Dirac Hamiltonian without disorder : Hq = ict • V + mas. 
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III. SADDLE POINT APPROXIMATION 



At E = the DOS of the pure Dirac fermions is always vanishing. Near m — ~ the point where the pure 
system closes the energy gap - the above supersymmetric effective action enables us to recover the average DOS of 
the random system. From the uniform solution for the SP equation (6S = 0) in the matrix argument Qr we obtain 
a semicircular DOS with radius TOc — 2e^'"l3 with respect to its dependence on the mass m jlj]. Thus this predicts, 
within a large- approach, a finite bandwidth for the states between two critical points m — ±mc. 

In order to work with localized states, a non-uniform solution for the SP equation for Qr must be sought, as in 
the approach by Cardy jTsj and Brezin jl^, who considered a random Schrodinger Hamiltonian using the bosonic 
replica trick. The situation is more complicated in the case of our supersymmetric field theory because of the presence 
of anticommuting variables. To avoid the difficulty of solving a differential equation with both commuting and 
anticommuting field components, only a schematic discussion in terms of the classical field equation (SP equation) is 
given in this section. From these qualitative arguments we then devise a procedure to obtain the large m behavior of 
the field theory. 



In order to obtain a non-uniform (soliton-like) solution for the SP equation of our action S, Eq. (2.4), it is convenient 
to work in the continuum limit of this model. For this purpose we introduce the characteristic length scale a/^^g on 
the J2r- Then the continuum limit, aj ^fg — > 0, can be taken. Our action thus transforms as follows 

5=y^(Tr2g? + ---) ^ J {Tv.Ql + ■ ■ ■)d\. (3.1) 

The SP method leads to the equation 

o-gTQrTCTs = g{Ho + iujcTo - 2tQt):^^. (3.2) 

We notice that we can replace Hq by ma^ -I- icr • V. Then the SP equation can be expanded in powers of 1/m because 
we are only interested in the large m behavior. To this end we pull out a factor ma^ 

cr3TQrTa3 = g— (1 + —a ■ Vers + i—as TQras)'^ (3.3) 

mm mm 

and multiply by us from both sides, which gives 

TQrT = ^(1 + — a • V(73 + z— as - — rQTa3)-V3. (3.4) 
mm mm 

Using the notation Qr = rQr-r and taking the limit uj — 0, the SP equation reads, in leading order in powers of 1/m 
and in lowest non-trivial order in VQ,- 



Qr ~ — = — — Qr H -^QrCTsQr (^ ' VCT3 QctsCT • V)rr + ' ' ' • (3.5) 

m mlm 



2 ^ 4 ^ ^ 2 



In the expansion we have taken a second order term in Q as well as the second order gradient term (the first order 
term being traceless). The constant term can be removed by introducing the shifted new field 

Qr — Qr — —o's 
m 

and by keeping only the leading order in 1 /m for each power of Qr- Then the SP equation can be written schematically 
(omitting the Pauli matrices) as 



Qr = — Qr + — (V'Q)r, (3.6) 



where a and h are numbers. This equation can be interpreted as the SP equation for the fonowing effective action 



Seff = j (^Q'"5^Q'-5^(VQ)^)^'-- (3-7) 
We are interested in finding the dependence of the effective action on the physical parameters m and g. For this 



purpose we use as an ansatz for the non-linear SP equation, Eq. (3.6) 
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9 V9 

where Q is now some dimensionless function. Moreover, we introduce the rescaled two-dimensional space variable 

n — —=r. 

V9 



Then we obtain from the differential equation (3.6) a non-linear differential equation for Q{R), where the coefficients 
are numbers 



9 

Furthermore, the effective action reads in terms of Q{R) and 



(3.8) 



m 
9 



Since the integral is now just a number (it depends, of course, on the function Q(i?) which needs to be determined 
as a solution of the differential equation (3.8)), the action is proportional to rri^/g. This implies that for large m the 
DOS must be of the form 

p(m) ^ poe-""*'/^. 

The parameter a and the coefficient po of the exponential factor will be evaluated below. 



IV. BEHAVIOR FOR A LARGE DIRAC MASS 



Having determined the functional dependence for the DOS via a saddle point approximation, in the following we 
shall use a more direct route by focusing on the large-m limit, or on the tail corrections, without actually using the 
saddle point approximation. The approach we use is based on an expansion of the action S in powers of 1/m up 
to lowest order in this parameter. It will result in an expansion in the parameter g/m? , so the following calculation 
holds good for the tails of the localized states in the DOS, and for weak disorder. 

For this purpose, we ca n se parate the off-diagonal contribution of i/o, namely ffp = iu • V, from the logarithmic 
term of the action of Eq. (|]J), as follows: 

Hq + iujOQ - 2tQt = Bq + H'a 
Ho + iujao + 2iTPT = Bp + H^. (4.1) 

Then we expand the logarithmic term of the action around the diagonal part of the Hamiltonian. We obtain, e.g., for 
the Q-dependent expression to leading order in the off-diagonal term 

lndet(i/o + it^CTQ - 2tQt) = lndct(mcr3 + icuao - 2tQt) + Tr(i/oS^^). (4.2) 

Thus, for the complete expression we get, from the expansion in terms of the off-diagonal part Hq 

e-^ = e-^' [1 - Tr(i7^B-i) + Ti{H[,B-^) + 4Tr(rerBp VOrS-i) + • • •] (4.3) 
with the unperturbed part of the action 

S' ^- ^(TrsQ^ + Ti^P^ + 2Ti2QrOr) - lndet(S,) + Indet(Bp). (4.4) 

9 J. 

This is a local action and all non-local expressions appear in the expansion terms. This means that the unperturbed 
"weight" e""^ factorizes on the lattice. In particular, we obtain for the diagonal part of the Green's function in leading 
order the expression 
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= J Qr„riro)VPr,VQr„Ver,VQro J] j V{r)VPrVQrVQrVQr, (4.5) 

r^ro 

where S" = i X;r(Tr2Qr + Tr2 if + 2X1-26^6^) and V{r) is the r-dependent part of the integrand (i.e. the ratio of the 

two determinants times the corresponding factor of e""^ ). Due to the supersymmetric representation and because 
terms involving Bp^ and Bq^ are local operator at site r, we have the following relation ||l5[| : 

V{r)VPrVQrVQr'DQr = 1. (4.6) 

Thus we are left with the computation of the integral at the site tq, namely: 

^ det2 (mcTa + iuGi^ + i2TProT) 

""^ det2(m(T3 + ibJGo — 2TQroT) 

1 

9 

For the Hermitean matrices Qro: Pro we use the following parametrizations 

90 + 93 qi-iq2\ p ^ ( PQ+ P3 Pi - ip2 
qi + iq2 go - 93 y ' \Pi+ iP2 Po - P3 

This gives, for instance 

Tt2QI^ + Tr2P2 ^ 2{pl +pl+pl+pl + ql + ql + ql + ql). 
We can immediately perform the P^o "integration (see Appendix A), obtaining 

j det2(ma3 + iuja^ + i2TProT) exp{-^{TT2Pro)T^Pro = -{^g fi^^ + + 2g). 



exph-(Tr2Q^,, + Tr2P4 + 2Tr2e.„e,.J]I?p.„I?Q.„De,„I?e.„. (4.7) 



(4.8) 



The integration over the Grassmannian fields O and Q contributes a factor {^)'^- Then the DOS reads, according to 
Eqs. (^ and (|2j|) 

p(m) = -4-^ lim (w^ + + 2g)Im [ — ^e-2('?o+9?+'?2+9i)/9dqp . . . ^^3. (4.9) 

TT-^gJ c^^o J det2[Bg) 

At this point, we have to carry out the integration with respect to Qro- This integral is evaluated in Appendix B and 
gives, finally 

p{m) = -^(1 + mV25)'e-'"'/45. (4.10) 



The typical form of the DOS, as a function of m and at the Fermi energy (i.e. £' = 0) is shown in Figure 1 (where a 
comparison with the corresponding extended-states (uniform SP solution) semi-circular DOS is also shown). We can 
now see that the DOS becomes, in the limit of large m and as anticipated, a function of m? / g which characterizes 
the tails of the localized electron states distribution. These tails become more pronounced as the disorder, and thus 
g, increases. The exponential dependence on m? / g is characteristic of a soliton-like solution of the SP approximation. 
In particular, we are now able to determine the parameters of the SP approximation of the previous Section as 

a = 1/4, Po ^ ^_(l + mV25)' 



without solving the SP equation (3^) directly. 
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V. CONCLUSIONS 



Dirac femiions have two bands with an energy gap proportional to the mass m. We evaluate the DOS at E — 0, 
i.e. in the middle of the gap, whilst varying the mass m. This means that we measure the DOS of the tails which 
develop inside the gap and find that the DOS acquires the asymptotic form (1 + rr? I2g)^e~'^ 1'^^ for small values 
of g/m?. Thus the behavior is controlled by the the exponential part e^™ Z**^. Conversely, in the case of ordinary 
Schrddinger particles subject to a Gaussian white-noise potential, Cardy Q and (with a different prefactor) Brezin 
and Parisi jlj] found tails in the energy E for the DOS which follow the more simple exponential decay. In our case, 
the Gaussian decay of the DOS with m reflects the Gaussian distribution of the Dirac mass. With other distributions 
for the disorder, we would expect also other forms for the DOS tails. 

We have addressed the question of the existence of localized states created by a random mass in the gap of the 
Dirac fermions' model in 2D. Our findings are in reasonable agreement with previous calculations for this model based 
on the homogenous SP method |^ and for 2D Dirac fermions with a random energy term 0. It is interesting that 
the expansion parameter is here g/w? (i.e. that the expansion is valid for weak randomness g), whereas localized 
states are usually related to strong randomness. Of course, it is crucial that also a relatively large mass m is present 
in order to obtain a valid expansion parameter. This regime is complementary to what we have in the homogeneous 
saddle point approximation where weak randomness is considered at any mass m. In the latter we have critical points 
TO = ±mc = ±26"^^/^. It would be interesting to approach these critical points, which are related to mobility edges 
with extended states for small |to|, using the non-local method developed in this paper. This could contribute, e.g., 
to the long-standing problem of the divergence of the localization length ^ (to) as the mobility edge is approached. 

Acknowledgements 

One of us (GJ) is grateful to Prof. K.B. Efetov and to the DFG via Sonderforschungsbereich 237 (Unordnung und grosse 
Fluktuationen, Projekt A8: " Fluktuationen in mesoskopischen und ungeordneten Systeme") at the Ruhr-Universitat Bochum, 
for partial support during the initial stage of this work. He is also grateful to Prof. P. Fulde and to the Max-Planck-Gesellschaft 
for support at the MPI fiir Physik Komplexer Systeme in Dresden. 

APPENDIX A 

In this Appendix we work out the P-block integration. With the expression 

^ y —2ip2 — 2pi lU) + 2ip2 — TO — 2ipQ J 

so that 

- det2Pp(ro) = {uj + 2ps^ro? + (m + 2ipo,ro)^ + 4pL„ + 4^? 

we see that the integration over the 2x2 P-block is easily carried out. Bearing in mind that Tr2(P^g) 
get 

- j dpidp2det2Ppe-f (^''+^'') = ^ ^ 3^3)' + (m + 2ipo)% 



(A.l) 
(A.2) 

(A.3) 



dp3[|<?' + \9{{io + 2p3)' + (m + 2ipof)]e 



+ -5(^2 ^ ^ 2ip^f )]{-g)-^ + 7r5(|)3 



(A.4) 



and finally 



.9' + |.9(^' + ("^ + 2iPo)')](|.g)^ + (|)^^^}e 



C^gnu^ + m^ + 2g]. 



(A.5) 
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APPENDIX B 

At this point we carry out the integration over the elements of the Q-block. We first of aU notice that 

7-. • r, ^ ( iui — Iq-i + m — 2on — 2ff2 — Vlq\ \ /n , \ 

i3„ = mcTs + iwcTn - 2rQT ^ • n (B.l) 



so that det2i3q = (iw — 2^3)^ — (m — 2go)^ + 4(7| + 4(7 J. Then we can write 

4 ^ + g| + g3 - (go -f?-^+ '^g3 

det2i?, (q2 + q2 + ^2 _ _ ™)2 _ <^)2 + ^2^2 ' 



(B.2) 



The contribution to the imaginary part of Eq. (4.9) comes from the integral 

T / ,222. 'L,2 .-^2 ^^^-^'^^+-+^^'^/^rfgo---rfg3. 

47 (g? + gi + gi - (go - f ) - "t) + ^ gi 

First of all, we perform the angular integration of the two-dimensional integral over qi and (72. This gives, with the 
definition q"^ = qf + g| 

^27r / / M , -e-^^^'+^o+^'^^/^dqdqodqs. (B.3) 

4 7 70 (g2+g|-(go-f)2-f)2+c.2g2 ^ ^ 

Next we evaluate the (/-integration for lu ^ 

2 r ge-^^^/^ 

'^'^'X ((z2+(7|-((Zo-m/2)2-c.W-f^2^| 

~ ||(?3|e^(^3-(.o-™/2)^-^V4)/.0((^^ _ „,/2)2 + ^2/4 _ ^2)^ 

and carry out the q^ integration 

|g3|e((go - m/2f + Lu^/A - ql)dq^ - ((jo - m/2f + u? jA. 

Finally, the (70-integration yields 

[((?o - m/2)2 + wV4]e-2(?o^+(*-m/2)=)/s^g^ ^ V^e-mV49(^ + ^2/2 + 20^"). 

16 



The combination all these results leads to an expression for Eq. (B.3) as 



- /4s(5-t-™V2 + 2c.2). (B.4) 
64 
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Figure Caption 

Density of states as a function of the average mass m at energy E — and disorder strength g = 2. The g/rn^ 
expansion (full curve) is compared with the result of the large- iV limit (homogeneous saddle point solution) of Ref. 
||12| (broken curve). 
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